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1. Origin of Higher Twist Terms 

The quark distributions in the proton and neu- 
tron are parametrized as Parton Distribution 
Functions (PDFs) obtained from global fits to 
various sets of data at very high energies. These 
fits are done within the theory of Quantum Chro- 
modynamics (QCD) in either leading order (LO) 
or next to leading order (NLO). The most impor- 
tant data come from deep- inelastic e/ fj, scattering 
experiments on hydrogen and deuterium, and 
and experiments on nuclear targets. In pre- 
vious publications |2I3I4| we have compared the 
predictions of the NLO MRSR2 PDFs to deep- 
inelastic e//z scattering data jjj on hydrogen and 
deuterium from SLAC, BCDMS and NMC. In or- 
der to get agreement with the lower energy SLAC 
data for Fi and R down to Q 2 = l GeV 2 , and at 
the highest values of x (x = 0.9), we found that 
the following modifications to the NLO MRSR2 
PDFs must be included. 

1. The relative normalizations between the 
various data sets and the BCDMS system- 
atic error shift must be included |2I3| . 

2. Deuteron binding corrections need to be ap- 
plied and the ratio of d/u at high x must be 
increased as discussed in ref. pj]- 

3. Kinematic higher- twist originating from 



target mass effects jS| are very large and 
must be included. 

4. Dynamical higher-twist corrections are 
smaller but also need to be included |2I3| . 

5. In addition, our analysis including QCD 
Next to NLO (NNLO) terms shows that 
most of the dynamical higher-twist correc- 
tions needed to fit the data within a NLO 
QCD analysis originate from the missing 
NNLO higher order terms. 

Our analysis shows that the NLO MRSR2 PDFs 
with target mass and NNLO higher order terms 
describe electron and muon scattering F2 and R 
data with a very small contribution from higher 
twists. Studies by other authors |xj also show 
that in NNLO analyses the dynamic higher twist 
corrections are very small. If (for Q 2 > 1 GeV 2 ) 
most of the higher-twist terms needed to obtain 
agreement with the low energy data actually orig- 
inate from target mass effects and missing NNLO 
terms (i.e. not from interactions with spectator 
quarks) then these terms should be the same in 
and e//i scattering. Therefore, low energy data 
should be described by the PDFs which are fit to 
high energy data and are modified to include tar- 
get mass and higher-twist corrections that fit low 
energy e//i scattering data. However, for Q 2 < 
1 GeV 2 additional non-perturbative effects from 
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Figure 1. Electron and muon F 2 data (SLAC, BCDMS, NMC, HI 94) used in our GRV98 £ w fit compared 
to the predictions of the unmodified GRV98 PDFs (LO, dashed line) and the modified GRV98 PDFs fits 
(LO+HT, solid line); [a] for F% proton, [b] for F 2 deuteron, and [c] for the Hf and NMC proton data at 
low x. 



spectator quarks must also be included 

2. Previous Results with GRV94 PDFs 

and x w 

In a previous communication .5] we used a 
modified scaling variable x w and fit for modifi- 
cations to the GRV94 leading order PDFs such 
that the PDFs describe both high energy and low 
energy e/[x data. In order to describe low energy 
data down to the photoproduction limit (Q 2 = 0) , 
and account for both target mass and higher twist 
effects, the following modifications of the GRV94 
LO PDFs are need: 

1. We increased the d/u ratio at high x as de- 



scribed in our previous analysis [2]- 

2. Instead of the scaling variable x we used the 
scaling variable x w = (Q 2 + B)/(2Mv + A) 
(or =x(Q 2 + B)/{Q 2 + Ax)). This mod- 
ification was used in early fits to SLAC 
data ^Oj- The parameter A provides for 
an approximate way to include both target 
mass and higher twist effects at high x, and 
the parameter B allows the fit to be used all 
the way down to the photoproduction limit 

(Q 2 =o). 

3. In addition as was done in earlier non-QCD 
based fits to low energy data, we mul- 
tiplied all PDFs by a factor K=Q 2 / (Q 2 
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+C). This was done in order for the fits to 
describe low Q 2 data in the photoproduc- 
tion limit, where F 2 is related to the photo- 
production cross section according to 

, , 4ir 2 a EM ^ 0.112m& GeV 2 
<*{1P) = n2 F * = 752 F2 



4. Finally, we froze the evolution of the GRV94 
PDFs at a value of Q 2 = 0.24 (for Q 2 < 
0.24), because GRV94 PDFs are only valid 
down to Q 2 = 0.23 GeV 2 . 

In our analyses, the measured structure func- 
tions were corrected for the BCDMS system- 
atic error shift and for the relative normaliza- 
tions between the SLAC, BCDMS and NMC 
data |2I3| . The deuterium data were corrected 
for nuclear binding effects |2I3| . A simultane- 
ous fit to both proton and deuteron SLAC, NMC 
and BCDMS data (with x > 0.07 only) yields 
A=1.735, B=0.624 and C=0.188 (GeV 2 ) with 
GRV94 LO PDFs (x 2 = 1351/958 DOF). Note 
that for x w the parameter A accounts for both 
target mass and higher twist effects. 

3. New Analysis with £ w , G D and GRV98 
PDFs 

In this publication we update our previous 
studies, [S] which were done with a new improved 
scaling variable £ w , and fit for modifications to 
the more modern GRV98 LO PDFs such that 
the PDFs describe both high energy and low 
energy electron/muon data. We now also in- 
clude NMC and HI 94 data at lower x. Here 
we freeze the evolution of the GRV98 PDFs at 
a value of Q 2 — 0.8 (for Q 2 < 0.8), because 
GRV98 PDFs are only valid down to Q 2 = 0.8 
GeV 2 . In addition, we use different photoproduc- 
tion limit multiplicative factors for valence and 
sea. Our proposed new scaling variable is based 
on the following derivation. Using energy mo- 
mentum conservation, it can be shown that the 
factional momentum £ = {p z + Po) / (P z + Po) car- 
ried by a quark of 4-mometum p in a proton tar- 
get of mass M and 4-momentum P is given by £ 
= xQ' 2 /[0.5Q 2 (l + [1 + (2Mx) 2 /Q 2 ] 1 / 2 )], where 



2Q' 2 = [Q 2 +Mf 2 -M l 2 ] + [{Q 2 +M f 2 -M l 2 ) 2 + 
4Q 2 (M> 2 + P 2 )] 1 ' 2 . 

Here Mi is the initial quark mass with average 
initial transverse momentum Pt and Mf is the 
mass of the quark in the final state. The above 
expression for £ was previously derived |H] for the 
case of Pt — 0. Assuming M, = we use instead: 

£ w = x{Q 2 + B + M/ 2 )/(0.5Q 2 (l + [1 + 
(2Mx) 2 /Q 2 ] 1 / 2 ) + Ax) 

Here Mf=0, except for charm-production pro- 
cesses in neutrino scattering for which Mf =1.5 
GeV. For £ w the parameter A is expected to be 
much smaller than for x w since now it only ac- 
counts for the higher order (dynamic higher twist) 
QCD terms in the form of an enhanced target 
mass term (the effects of the proton target mass 
are already taken into account using the exact 
form in the denominator of £ w ). The parame- 
ter B accounts for the initial state quark trans- 
verse momentum and final state quark ef fective 
AMf 2 (originating from multi-gluon emission by 
quarks). 

Using closure considerations (e.g. the 

Gottfried sum rule) it can be shown that, at 
low Q 2 , the scaling prediction for the valence 
quark part of F 2 should be multiplied by the 
factor K=[1-G 2 D (Q 2 )][1+M(Q 2 )] where G D = 
l/(l+<3 2 /0.71) 2 is the proton elastic form fac- 
tor, and M(Q 2 ) is related to the magnetic elastic 
form factors of the proton and neutron. At 
low Q 2 , [l-G 2 D {Q 2 )] is approximately Q 2 /(Q 2 
+C) with C = 0.71/4 = 0.178 (versus our fit 
value C=0.18 with GRV94). In order to satisfy 
the Adler Sum rule Q] we add the function 
M(Q 2 ) to account for terms from the mag- 
netic and axial elastic form factors of the nu- 
cleoli). Therefore, we try a more general form 
K va i ence =[l-G 2 D (Q 2 )}[Q 2 +C 2v ]/[Q 2 +Ci„], and 
K sea =Q 2 /(Q 2 +C sea ). Using this form with the 
GRV98 PDFs (and now also including the very 
low x NMC and HI 94 data in the fit) we find 
A=0.419, 5=0.223, and Ci„=0.544, C 2u =0.431, 
and C sea =0.380 (all in GeV 2 , X 2 = 1235/1200 
DOF). As expected, A and B are now smaller 
with respect to our previous fits with GRV94 
and x w . With these modifications, the GRV98 
PDFs must also be multiplied by iV=1.011 to 
normalize to the SLAC F 2p data. The fit (Fig- 
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Figure 2. Comparisons to proton and iron data not included in our GRV98 £ w fit. (a) Comparison of 
SLAC and JLab (electron) F 2p data in the resonance region (or fits to these data) and the predictions 
of the GRV98 PDFs with (LO+HT, solid) and without (LO, dashed) our modifications, (b) Comparison 
of photoproduction data on protons to predictions using our modified GRV98 PDFs. (c) Comparison of 
representative CCFR and F M charged-current differential cross sections |4I14| on iron at 55 GeV and 
the predictions of the GRV98 PDFs with (LO+HT, solid) and without (LO, dashed) our modifications. 



ure ^) yields the following normalizations rel- 
ative to the SLAC F 2p data {SLAC D =0.986, 
BCDMS P =0.964, BCDMS D =0.9M, 
NMC P =1.00, NMC D =0.993, m P =0.977, and 
BCDMS systematic error shift of 1.7). (Note, 
since the GRV98 PDFs do not include the charm 
sea, for Q 2 > 0.8 GeV 2 we also include charm 
production using the photon-gluon fusion model 
in order to fit the very high v HERA data. This 
is not needed for any of the low energy compar- 
isons but is only needed to describe the highest 
v HERA electro and photoproduction data). 

Comparisons of predictions using these modi- 
fied GRV98 PDFs to other data which were not 
included in the fit is shown in Figures |21 and [3] 
From duality 15 considerations, with the £ w scal- 



ing variable, the modified GRV98 PDFs should 
also provide a reasonable description of the av- 
erage value of F% in the resonance region. Fig- 
ures EJa) and |3{ a) show a comparison between 
resonance data (from SLAC and Jefferson Lab, 
or parametrizations of these data )16j ) on pro- 
tons and deuterons versus the predictions with 
the standard GRV98 PDFs (LO) and with our 
modified GRV98 PDFs (LO+HT). The modified 
GRVB98 PDFs are in good agreement with SLAC 
and JLab resonance data down to Q 2 = 0.07 (al- 
though resonance data were not included in our 
fits). There is also very good agreement of the 
predictions of our modified GRV98 in the Q 2 = 
limit with photoproduction data on protons and 
deuterons as shown in Figure CJb) and^b). In 
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Figure 3. Comparisons to data on deutrerium which were not included in our GRV98 £ w fit. (a) Compari- 
son of SLAC and JLab (electron) F%d data in the resonance region and the predictions of the GRV98 PDFs 
with (LO+HT, solid) and without (LO, dashed) our modifications, (b) Comparison of photoproduction 
data on deuterium to predictions using our modified GRV98 PDFs (including shadowing corrections), 
(c) The shadowing corrections that were applied to the PDFs for predicting the photoproduction cross 
section on deuterium. 



predicting the photoproduction cross sections on 
deuterium, we have applied shadowing correc- 
tions 2Q; as shown in Figure Etc)- We also com- 
pare the predictions with our modified GRV98 
PDFs (LO+HT) to a few representative high en- 
ergy CCFR and 77^ charged-current differential 
cross sections |4I14) on iron (neutrino data were 
not included in our fit). In this comparison we 
use the PDFs to obtain F 2 and xF 3 and correct 
for nuclear effects in iron jS] . The structure func- 
tion 2x F\ is obtained by using the R WO rid fit from 
reference [5] . There is very good agreement of our 
predictions with these neutrino data on iron. 

In order to have a full description of all charged 
current and processes, the contribution 
from quasielastic scattering ^7] must be added 
separately at x = 1. The best prescription is to 
use our model in the region above the first reso- 



nance (above Vy=1.35 GeV) and add the contri- 
butions from quasielastic and first resonance |18j 
(VF=I.23 GeV) separately. This is because the 
W = M and W—\.2i GeV regions are dom- 
inated by one and two isospin states, and the 
amplitudes for neutrino versus electron scatter- 
ing are related via Clebsch-Gordon rules JHj in- 
stead of quark charges (also the V and A cou- 
plings are not equal at low W and Q 2 ). In the 
region of higher mass resonances (e.g. W=1.7 
GeV) there is a significant contribution from the 
deep-inelastic continuum which is not well mod- 
eled by the existing fits |IB| to neutrino resonance 
data (and using our modified PDFs should be bet- 
ter). For nuclear targets, nuclear corrections 8 
must also be applied. Recent results from Jlab 
indicate that the Fe/D ratio in the resonance re- 
gion is the same as the Fe/D ratio from DIS data 
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for the same value of £ (or £ w ). The effects of 
terms proportional to the muon mass and F4 and 
F 5 structure functions in neutrino scattering are 
small and are discussed in Ref . |17I19| . In the fu- 
ture, we plan to investigate the effects of including 
the initial state quark Pp in £ w , and institute fur- 
ther improvements such as allowing for different 
higher twist parameters for u, d, s, c, b quarks 
in the sea, and the small difference (expected in 
the Adler sum rule) in the K factors for axial and 
vector terms in neutrino scattering. In addition, 
we can multiply the PDFs by a modulating func- 
tion |10II2| A(W,Q 2 ) to improve modeling in the 
resonance region (for hydrogen) by including (in- 
stead of predicting) the resonance data JS| m the 
fit. We can also include resonance data on deu- 
terium |16| and heavier nuclear targets in the fit, 
and low energy neutrino data. Note that because 
of the effects of experimental resolution and Fermi 
motion |2] (for nuclear targets), a description of 
the average cross section in the resonance region 
is sufficient for most neutrino experiments. 

The current analysis assumes that the ax- 
ial and vector structure functions are equal at 
all Q 2 . However, at very low Q 2 , the vector 
structure function must go to zero, while the 
axial-vector part is finite. We are currently (Au- 
gust 2003) in the process of including low energy 
data from Chorus (on Pb) in our fit, in order to 
constrain the low Q 2 axial-vector contribution. 
As for the vector case, the form of the fits is 
motivated by the Adler sum rule for the axial- 
vector contribution as follows: K va i ence - ax = 
[1-F% (Q 2 )] [Q 2 +D 2v - ax ] I [Q 2 ( +D lv _ ax ], 

and K sea — ax — ~^~-&sea — ax)/{.Q ~\~ $ sea— ax ) • 

Here F A = -1.267/(1+Q 2 /1.00) 2 . 



4. Appendix 

In leading order QCD (e.g. GRV98 LO PDFs), 
F2 for the scattering of electrons and muons 
on proton (or neutron) targets is given by the 
sum of quark and anti-quark distributions (each 
weighted the square of the quark charges) : 



2xF 1 (x) 



2 r xqi(x) + xqt(x)] 
F 2 {x)(l+4Mx 2 /Q 2 )/(l 



(1) 



Here, R w (x,Q 2 ) is parameterized [5] by: 
0.0635 

= lo 9 (Q 2 /0M) 9{x > Q) 
0.5747 0.3534 



Q 4 + 0.09 



(3) 



where 6 = 1. 



12Q 2 

Q' 2 +i.o 



0.125 2 



0.125^+2^' 

The R w function provides a good description 
of the world's data in the Q 2 > 0.5 and x > 0.05 
region. Note that the R w function breaks down 
below Q 2 = 0.3. Therefore, we freeze the function 
at Q 2 = 0.35 and introduce the following function 
for R in the Q 2 < 0.35 region. The new func- 
tion provides a smooth transition from Q 2 = 0.35 
down to Q 2 = by forcing R to approach zero at 
Q 2 = as expected in the photoproduction limit 
(while keeping a 1/Q 2 behavior at large Q 2 and 
matching to R w at Q 2 = 0.35). 



R = 3.207 x 



Q 2 



Q 4 



1 



R w {x,Q 2 =0.35). (4) 



In the comparison with CCFR charged-current 
differential cross section on iron, a nuclear correc- 
tion for iron targets is applied. We use the fol- 
lowing parameterized function, f(x) (fit to exper- 
imental electron and muon scattering data for the 
ratio of iron to deuterium cross sections), to con- 
vert deuterium structure functions to (isoscalar) 
iron structure functions 



f(x) = 1.096- 0.364a; - 0.278e 
+2.772a; 14 - 417 



-21.94a: 



(5) 



For the ratio of deuterium cross sections to 
cross sections on free nucleons we use the follow- 
ing function obtained from a fit to SLAC data 
on the nuclear dependence of electron scattering 
cross sections 0]. 

/ = (0.985 ± 0.0013) x (1 + 0.422a; - 2.745a; 2 
+7.570a; 3 - 10.335a; 4 + 5.422a; 5 ). (6) 

This correction is only valid in the 0.05 < x < 
0.75 region. In neutrino scattering, we use the 
same nuclear correction factor for F 2 , xF% and 
2xF t . 

The d/u correction for the GRV98 LO PDFs is 
obtained from the NMC data for F^/F 2 P . Here, 
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Eq. is used to remove nuclear binding effects 
in the NMC deuterium F 2 data. The correction 
term, 8{d/u)(x) is obtained by keeping the total 
valence and sea quarks the same. 



5{d/u) 



-0.008f7 + 0.0506a; + 0.0798a; 2 , (7) 



where the corrected d/u ratio is (d/u)' = (d/u) + 
8{d/u). Thus, the modified u and d valence dis- 
tributions are given by 



!/.,, 



d' 



d v + u v S(d/u) 



1 



(8) 
(9) 



The same formalism is applied to the modified u 
and d sea distributions. Accidently, the modified 
u and d sea distributions (based on NMC data) 
agree with the NUSEA data in the range of x be- 
tween 0.1 and 0.4. Thus, we find that any futhcr 
correction on sea quarks is not necessary. 
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